EPJ manuscript No. 

(will be inserted by the editor) 



On the surface tension of fluctuating quasi-spherical vesicles 



o 



o 



C. Barbetta^'', A. Imparato^, and J.-B. Fournier^ 

^ Laboratoire Matiere et Systemes Complexes (MSG), UMR 7057 CNRS & Universite Paris Diderot-Paris 7, 

Bat. Condorcet, CC 7056, 75205 Paris, France. 
^ Department of Physics and Astronomy, University of Aarhus, Ny Munkegade, Building 1520, 

DK-8000 Aarhus C, Denmark 

Received: date / Revised version: date 

Abstract. We calculate the stress tensor for a quasi-spherical vesicle and we thermally average it in order 
to obtain the actual, mechanical, surface tension r of the vesicle. Both closed and poked vesicles are 
considered. We recover our results for r by differentiating the free-energy with respect to the proper 
projected area. We show that r may become negative well before the transition to oblate shapes and that 
it may reach quite large negative values in the case of small vesicles. This implies that spherical vesicles 
may have an inner pressure lower than the outer one. 

PACS. 87.10.-e Biological and medical physics: general theory and mathematical aspects - 87.16. dj Mem- 
branes, bilayers, and vesicles: dynamics and fluctuations 
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1 Introduction 

Biological membranes are nanometer-thick sheets com- 
posed essentially of a bilayer of phospholipids [1] . Although 
very thin, they form a barrier for ions and larger molecules. 
They are thus essential for cells, forming a protection from 
the external media and preventing diffusion. Moreover, 
they are fluid and very flexible, exhibiting a strong resis- 
tance to area change but no shear resistance. As a conse- 
quence, they present large thermal fluctuations [5]. 

Bilayers composed of one single type of phospholipid 
are usually used to study the physical properties of mem- 
branes. For this purpose, unilamellar vesicles, consisting 
of one closed bilayer, are widely used in experiments. Dif- 
ferent techniques allow to obtain these vesicles, whose size 
range from a few tens of nanometers to a hundred of mi- 
crometers (giant vesicles) [3l|4j[5]. They appear in several 
fluctuating shapes, depending on the enclosed volume V, 
on the total area A and on the area difference between the 
monolayers |6l|7j. 

For liquid interfaces, the surface tension is a constant 
depending only on the nature of the phases in contact. 
The situation for membranes is quite different. Membranes 
consist of a fixed number of insoluble lipids possessing an 
equilibrium density. Mechanically, the tension depends on 
the lipid density but also on the curvature of the mem- 
brane [51151110]. Nonetheless, for fluid membranes, the con- 
cept of surface tension is threefold: (i) Experimentally, 
one measures a tension r from the quadratic dependence 
in the wavenumber of the fluctuation spectrum [5]. (ii) 
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Theoretically, one can introduce in the Hamiltonian a La- 
grange multiplier a multiplying the total membrane area. 
In principle r can be deduced from a by renormalization 
procedures [mil2j . (More accurately, one may fix the total 
membrane area by a constraint |13j or, even better, keep 
the lipid density as a parameter [5].) (iii) Experimentally, 
one measures another tension, r, through the pressure dif- 
ference across the membrane by means of the Laplace rela- 
tion P^ITB] . In principle, this is the actual (macroscopic) 
tension, since it corresponds to the average force that is 
exerted tangentially to the average membrane's surface. 
In principle, r can be calculated by differentiating the to- 
tal free-energy with respect to the projected area [H] , but 
this is a tricky matter |16] . The confusion associated with 
these different tensions inspired many articles [Tn[T2l[T6l 

nziiiiiiiiiiniiii]. 

Recently, a new method was introduced to calculate 
T, based on thermally averaging the stress tensor of the 
membrane [TB] . It has the advantage to be directly related 
to the very definition of r. This method was applied to 
planar membranes and it was shown that the expression 
of T obtained in this way coincides with that obtained 
from carefully differentiating the total free-energy fl6| . 

The aim of this paper is to calculate the surface tension 
T of quasi-spherical vesicles by using the same method 
as in |16j . i.e., by averaging the stress tensor. There are 
several interesting questions that we shall address: 

1. What is the difference between r for a vesicle and r 
for a planar membrane? Is there a characteristic radius 
over which they coincide? 

2. How does the volume constraint influence the expres- 
sion for T? 
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Fig. 1. Parametrization in spherical coordinates of a vesicle 
(bold line) fluctuating around a reference sphere (dashed line). 
The inset shows the force exchanged through a cut that sepa- 
rates region 1 and region 2. 



3. Can one obtain r by differentiating the total free-energy 
with respect to the projected area? If yes, what does 
projected area mean in the case of a vesicle? 

4. Can T become negative (in this case the Laplace pres- 
sure would be inverted with respect to the normal sit- 
uation)? 

Our paper is organized as follows. In sec. [5] we pre- 
cisely define the system under study and we work out the 
membrane stress tensor in spherical coordinates. Then we 
calculate the various correlation functions involved in the 
vesicle's shape fluctuations. In sec. 121 we use these cor- 
relation functions to calculate r by thermally averaging 
the stress tensor. We discuss the relation between r and 
the membrane excess area. In sec. SI we succeed in recov- 
ering T by differentiating the free-energy with respect to 
the vesicle's and we discuss the influence of the vol- 
ume constraint. In scc.[5]wc summarize our results and we 
discuss the issue of negative tensions. 

2 Stress tensor and fluctuation spectrum 

2.1 Effective Hamiltonian and parametrization 

We consider a quasi-spherical vesicle with flxcd area A 
and fixed volume V . Its shape is parametrized in spherical 
coordinates by (see Fig. [T|) 

r = R[l + u{e,(t))]er, (1) 

where u <C 1 . A convenient choice for the reference sphere 
is i? = (|y/7r)^/"^, corresponding to the sphere of volume 
y|13j. The volume and area constraints read: 

V^-R^i del d4)[l + u{e,4))f sine , (2) 
^ JO Jo 



A = de d(l){l + u) 

Jo Jo 

X ^J{1 + uf sir? 6 + + sin^ 9 . (3) 

Here, and in the following, u,; = du/di, Uij = d^u/didj, 
where i, j G {9, 0}. Latin indices will denote either 9 or 
(j), not r. 

There arc three variants of the free-energy describ- 
ing the elasticity of a fluid vesicle: i) the "spontaneous 
curvature", or Helfrich model [12], ii) the "bilayer cou- 
ple" model [2511^ , and iii) the "area-difference elasticity" 
model [7l [25l[26] . the latter being the most accurate [7]. 
As shown by Seifert [T3], the three models are all equiv- 
alent in the quasi-spherical limit to the minimal Helfrich 
model, i.e., to the spontaneous curvature model with van- 
ishing spontaneous curvature [22] . We therefore adopt the 
latter, which corresponds to the free-energy, or effective 
Hamiltonian: H = 2k J dS , supplemented by the area 
and volume constraints. Here the integral runs over the 
vesicle's surface, C is the local mean curvature of the mem- 
brane and K is the bending rigidity constant. While the 
volume constraint is quite easy to implement, it is diffi- 
cult to handle the surface constraint exactly [13] . We shall 
therefore use the traditional approach, in which the area 
constraint is approximatively taken into account by means 
of a Lagrange multiplier a playing the role of an effective 
surface tension. As discussed in Ref. [13], the latter ap- 
proach gives correct results in the small excess area limit 
(in which we shall place ourselves) . 

Thus, the effective Hamiltonian we shall use is 

JdS{2KC^+a), (4) 

together with the volume constraint ([2]). Explicitly, in 
terms of u{9, (p), the Hamiltonian reads 

H = J d9 d4>h{u,{ui},{uij}) , (5) 

with [MIUS]: 

h ^ {2k + R^ a) sine 
+ 2 sin 9 [R^<7 u — k {u^^ csc^ 9 + ug cot 9 + ugg)] 

+ ]^sin9 [2i? V V? + {2k + R^(j) {ul + it^ csc^ 9) 

+ K [ug cot 9 + U^^ CSC^ 9)"^ + K Ugg{ugg + Au) 

+ 2K{ugg + 2u) {ug cot 9 + u^^ csc^ 9)] + 0{u^) . (6) 
2.2 Stress tensor components 

Let us consider an infinitesimal cut at constant longitude 
{(j) constant) separating a region 1 from a region 2 (see 
Fig. [T]). The normal to the projection of this cut onto the 
reference sphere is m = e^. By definition, the projected 
stress tensor S in spherical geometry relates linearly the 



C. Barbetta et al.: On the surface tension of fluctuating quasi-spherical vesicles 



3 



force dF that region 1 exerts on region 2 to the angu- 
lar length ds = dO of the projection of the cut onto the 
reference sphere: 

dF = S - mds = {Eg^ eg + + Sr^, e^) d9. (7) 

Likewise, for a cut at constant latitude {9 constant), with 
m = eg and ds = d(j), we have 

dF = S - mds ^ [Sgg eg + S^g 60 + Erg er) d(j). (8) 

For an oblique cut, dF is obtained by decomposing m 
along 60 and eg. 

The projected stress tensor associated with the Hamil- 
tonian of the minimal Helfrich model has been calculated 
in planar geometry |10j and in cylindrical geometry |27j 
(see also Ref. [5] for the original covariant formulation). 
The derivation in spherical geometry is lengthier but it 
follows the same route, which we sketch now. We consider 
a patch of membrane delimited by a closed curve, corre- 
sponding to a domain Q on the reference sphere enclosed 
by dn. The membrane within the patch is assumed to be 
deformed, at equilibrium, by means of a distribution of 
surface and boundary forces (and a distribution of bound- 
ary torques). To each point of the patch, we impose an 
arbitrary displacement 5a = Sur 6^ -I- Sag eg + 5a^ e^ that 
keeps constant along the boundary the orientation of the 
membrane's normal n (in this way the torques will pro- 
duce no work). On the one hand, the boundary energy 
variation reads (after integration by parts): 



5H = 



/ ds TTli 


hSi+i^ 


'dS2 





' dh 



dh 



dh 



duij ) duij 



.(9) 

where ds is the arc-length in the {9, </>) space, h is given 
by eq. m is the normal to dQ^ and Si G {i56',(5(/>} 
corresponds to the variation of dfi. On the other hand, 
the work of the force exerted through the boundary reads 



SH = ds Sa ■ S ■ m . 

Jon 



(10) 



By identifying eqs. (|9]) and ((T0| . one can obtain S. Details 
of this calculation, as well as all the components of the 
stress tensor are given in appendix |21 In particular, we 
obtain 

Sgg = ■^^R^asu\e{2 + 2u + ulcsc^9 -ufj 
+ K [u00 csc^ 9 - Ugg sin 9 

+ 2ug {ug^^ CSC 9 + Uggg siu 9)] 

+ K CSC 9 [2m0 — cos^ 9 — 2-U00 {1 + ug cot 9)] 
+ 4:Ku CSC 9 — Ugg sin^ 9) 

+ 2k [ugg (sin 9 + ug COS 9) + (2u — l)ug COS 9] I 

+ 0{u'). (11) 

2.3 Thermal averages and fluctuation spectrum 

In order to calculate the vesicle's mechanical tension r, 
we need to determine the thermal average (Egg). To this 



aim, we perform the standard decomposition of u(9, (j)) in 
spherical harmonics [T3l[28l[29] : 



"0,0 



with 



/Air 

-1)"<,„ and 



(12) 



E-EE 



(13) 



1=2 m = -l 



L being a high wavevector cutoff (see below). Note that 
the modes 1 = 1, which correspond to simple translations, 
are discarded. 

In terms of and up to order it^, eq. ^ takes the 
form: 



V = R^ 



An 



"0,0 



E 1"'^'^ 



(14) 



The volume constraint V = ^ttR^ (recall the definition of 
R) implies therefore [T^ : 



uo,o 



^E 

47r ^ 



With the help of the relation: 

Q2Y/n ^ g2y-m 



cot 9- ' 



-csc^6l- 



d9 d(j)^ d9^ 

and using eq. ([T5|). one obtains then |13l : 



(15) 



-i{i + \)Yr, (16) 



where 



Here, 



Hi = + 2) {P + 1 + a) 



k/R^ 



(17) 

(18) 
(19) 



is the reduced tension. Note that we have discarded in H 
a constant energy term, Sttk. 

We emphasize that negative values of a are allowed [13] . 
Indeed, the minimum of the Hamiltonian (|17p corresponds 
for a > —6 to a perfectly spherical vesicle (u/^m = 0, 
VZ > 2). The mean-field transition to an oblate shape oc- 
curs thus at fj = —6 (non harmonic terms being then 
needed to stabilize the system). 

Standard statistical mechanics yields then (M;.m) = 0, 
\Jl ^ 0, and 

{Ul,m ui,,n') = (-1)™^ <5u' <5™,-™' , (20) 
Hi 

where k^T = 1//3 is the temperature in energy units. 
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We may now calculate the fluctuation amplitudes. Us- 
ing eq. (j20[) and the addition theorem for spherical har- 
monics: 



m— — l 



21 + 1 



(21) 



we obtain 



47r 



47J- 



kBT ^ 21 



1=2 



H, 



(22) 



E 



fcRT ^ 2? 



4tt 



sm 



;=2 



4n ^ 

1=2 



+ {21 + 1) 
2i/; 



47r ^ 

;=2 



^(/ + 1)(2/ + 1) 
2i?/ 



etc. 



(23) 
(24) 
(25) 



The correlations of the other derivatives of u are given 
in appendix [BJ Note the interesting relation (u) = —{v?) 
(valid at second order in w), showing how the temperature- 
dependent fluctuations affect the mean-shape. 



2.4 Cutoff 

The large wavenumber cutoff L should be related with 
the smallest wavevector allowed, A w a~^, where a is a 
length comparable to the membrane thickness (i.e., it/ A 
of the order of a few times a). With spherical harmonics, 
however, this is not easy to implement. The requirement 
that we should recover the planar limit for large values of 
R will guide us. 

For a square patch of fluctuating membrane with refer- 
ence area Ap and periodic boundary conditions, the wave- 
vectors are quantified according to <7 = 2T:yJ Ap [ux, riy), 
where and Uy are integers and \q\ < A. The number of 
modes is then tt/I^/ {2tt/ y^A^)^ and the number of modes 
per unit area is 

AT /1 2 



An 



47r 



For the vesicle, we have 



modes 



47ri?2 



E (2^ + 1) = 



1=2 



(£-1) (L + 3) 
47ri?2 



(27) 



Asking that the number of degree of freedom per unit 
area (per lipid, in some sense) be the same in both case, 
we require these two quantities to be equal. Hence, we get 



{L - 1) (L + 3) = A^R^ 



(28) 



which gives L = \yA + R?A^ ~ Ij ( L^J is the integer part 
of x). In the limit R ^ A"^ , this gives simply L ~ AR. 



2.5 Validity of the Gaussian approximation 

Since our calculations are limited to ©(u^), wc should 
check, in principle, that higher order terms are negligi- 
ble. In practice this is not feasible. To check the smallness 
of u (which is especially critical in the case cr < 0) we shall 
require: 



2/ + 1 



47rK ^ (; - 1) (/ + 2) + 1 + a) 



To be definite, in the following we shall take 

U niax = 5 % . 



(29) 



(30) 



Note that the presence of the factor {P + I + a) in the 
denominator of eq. ((29)) . together with the condition / > 2, 
implies a G [—6,00], as already discussed. 

Solving condition (P^)) for the typical values A^^ ~ 
5nm, K = 25 k^T, and taking t/max = 0.05, we find 



cr > cr,. 



(31) 



almost independently of R. It follows that negative ten- 
sions cr are in fact within the validity range of our Gaussian 
approximation. 

An interesting control parameter is the excess area, 
defined by 



(A) Ap 



kjiT 



E' 



2/ + 1 



inn ^ P + I + a 



(32) 



Our validity condition pip yields a < Q!max, with Omax 
shown in fig. [2l One can see that Ofmax ~ ci -|- C2lni?, 
where ci and C2 are constants. Indeed, the sum in eq. ([5^ 
is dominated by the modes I = 2 and I = 3, the rest being 
well approximated for a = 0(1) by an integral propor- 
tional to In R. 



3 Actual mechanical surface tension 

Our aim is to calculate the mechanical tension r of the 
vesicle. Imagine replacing the fluctuating vesicle by a shell 
coinciding with its average shape. The actual mechanical 
tension r is the average force per unit length that is ex- 
erted tangentially to the shell's surface. Because of the 
spherical symmetry, we may evaluate it at any point in 
any direction. Let us consider at {9, (f) a cut with 9 con- 
stant of extension d(f). The component along eg of the force 
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2.5 I- 




i- i?[nml 

5C 

Fig. 2. Maximum excess area corresponding to \/Ju'^) < 0.05, 
which we take as the validity criterion of our Gaussian approx- 
imation. In abscissa is the vesicle's radius. Here, ~ 5nm 
and K = 25kBT. 



exchanged through the cut is on average {IJggd(j)). The 
length of the cut is on average + u) sinfld^). Hence, 



(So 



Rsme{l + (u)) 



(33) 



Since (u) = — (u^), because of the volume constraint, and 
since See = aRs'mO + 0{u), we obtain equivalently 



1 



i?sin6' 



{See) + a{u^) + 0{u^) . 



(34) 



Using eq. (fTTj) and the results obtained in sec. 12.31 and 
appendix [B] we obtain 



1=2 
L 



a 



E 



1) /(/ + !)(/ + 2) (2/ + 1) 
Hi 

1(1 + 1) (21 + 1) 



87ri?2 ^ P + l + a 



(35) 
(36) 



As expected this result is independent of point {6, (/>) where 
the calculation was made. 

We show in fig. |3] the behavior of r as a function of 
the excess area a, obtained numerically from cqs. (|36p 
and (j32p . This presentation is the most physical, because 
T and a are both measurable, while a is not. There are 
several salient points: 

1. The results for t deviate from the flat limit (i? ^ oo) 
essentially for R < 1 jiui. 

2. Negative and quite large values of r are indeed acces- 
sible. 

3. There exists a well defined excess area corresponding 
to T = 0. Since, as we shall show, the curves in fig. ^ 
are similar without the volume constraint, this value 
corresponds to the spontaneous excess area taken up 
by a poked vesicle (vanishing Laplace pressure). 

4. There is a plateau at large values of a, which probably 
corresponds to the actual transition to oblate shapes. 




a{%) 



Fig. 3. Mechanical tension r as a function of the excess area 
a, for AC = 10"^^ J, ksT = 4 X 10"^^ J and yl = (5nm)"\ The 
dashed lines correspond from the left to the right to vesicles of 
radiuses: R = 50 nm, R — 0.5 fim and R = 5 /im. The solid gray 
line corresponds to a flat membrane, eq. (|39}. The end-points 
indicate the limit beyond which our Gaussian approximation 
is no longer valid according to sec. 12.51 The curves for r given 
by eq. (49) (poked vesicles) are qualitatively similar. 



It occurs for f = tR'^/k < —6, i.e., below the mean- 
field threshold (see the discussion after eq. ([H])). The 
high symmetry phase (spherical vesicle) is thus stabi- 
lized by its entropic fluctuations, as one might have 
expected. 

We have also calculated the normal and orthogonal 
components of the tension. They vanish: {S^e) = and 
{S^e) = 0. While the latter result is obvious on symme- 
try grounds, the former one is interesting, implying that 
the shell mentioned above can indeed be considered as a 
purely tense surface. (This would probably not hold for 
a vesicle with a non-spherical average shape.) As a con- 
sequence, the Laplace law can be used without curvature 
corrections for a fluctuating quasi-spherical vesicle — if one 
uses T instead of a. Indeed, this could be expected from 
renormalisation arguments, since the Laplace law is exact 
(despite the curvature energy) for a perfectly spherical 
membrane |15| . 

It is interesting to examine r in the limit of large vesi- 
cles. In this case, the sum in eq. (j36[) may be replaced by 
an integral: 



87ri?2 
ksTA^ 
8^ 



^^^^ 1(1 + 1) (21 + 1) 



a 



P + l + a 
In I 1 



(37) 



O 



(38) 
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R !nml 



Fig. 4. Spontaneous excess area, corresponding to t = 0, as 
a function of the vesicle radius. The dashed line gives the flat 
membrane limit ao- The material parameters are the same as 
in flg. El 



The dominant term in eq. (|38l) correctly matches the ex- 
pression of T — (7 for flat membranes [16j . In this case, 
the relation between a and a is analytical |30| . yielding 
a = Kyl^/[exp(87r/3KQ;) — 1] and 



Tflat(a) 



kA^ (1 + a) kjiTA^ 



8n 



(39) 



which is shown in fig. [31 

Finally, we show in fig. |4] the excess area corresponding 
to a vanishing lateral tension (t = 0). Its value is very 
much radius dependent for R < 1/xm, but one recovers for 
R>2 /im the flat membrane limit [16] : 



ao 



In (87r/3K) 
87r/3K 



4 Derivation using the free-energy 

For a flat membrane, one may also obtain r by differ- 
entiating the free-energy with respect to the projected 
area Ap [T ^ I ^ I ^ ITC]. but there are two pitfalls HI]. One 
must: i) take the thermodynamic limit {Ap — >■ oo) only af- 
ter the differentiation, and ii) introduce a variation of the 
cutoff in order that the total number of modes remains 
constant during the differentiation, as discussed in [21j . 

Let us investigate the free-energy method in the case 
of quasi- spherical vesicles. The free-energy, F, is given by 



(41) 



the integral running over all the configurations of the vesi- 
cle. At the Gaussian level, H is given in terms of spherical 
harmonics by eq. (fTTj). and since r = Rer + Ru{9, (j)) e^, 
we may write (in agreement with Ref. |13|): 



V\r] = 



n( n RduiA ( n ^^^J , 

1=2 \m=0 / \m=l I 



(42) 



where the superscripts 7?. and I signify real part and imag- 
inary part, respectively. This measure corresponds to the 



so-called normal gauge, which is known to be correct for 
small fluctuations [13] . We note that the radius R of the 
reference sphere appears explicitly and that for each value 
of I, only half of the allowed values of m have to be con- 
sidered, as r is real. 

Performing the Gaussian integrals yields 



F = 47rR^a + kBTY^ 

1=2 



21 + 1 



In 



[l3Hi/R' 



(43) 



In order to obtain r, we must differentiate F with respect 
to the vesicle's "projected area" Ap. Which one, however? 
The area AttR^ of the reference sphere (i.e., the sphere 
having the same volume as the vesicle's), or the area of the 
vesicle's average shape, defined as A„i = 47r(i?(l + it))^? 
It will turn out that the former choice is the correct one. 
In a sense, this is natural because it corresponds to our 
parametrization. However, it is not that obvious, because 
the definition of r in eq. ([55]) involves the area of the 
average vesicle's shape. 

Let us thus pick Ap = AttR^. It is worth noticing that 

Hi depends on Ap only through a = aR^/n, yielding 



dHi 



With the choice: 



= il-l){l + 2)a. 



dF 1 dF 



dAp At: dR'^ 



(40) we obtain 



8ttR^ 



E 

1=2 



I {I + 1) (21 + 1) 
P + l + a 



(44) 



(45) 



(46) 



which is identical to the result obtained from the stress 
tensor approach, eq. (|36p . How about the pitfalls men- 
tioned above? First, we didn't take the thermodynamic 
limit before differentiating. (Actually, this would not be 
problematic, since the quantification of the modes does 
not involve the size of the system, like it is the case for pla- 
nar membranes.) Second, we kept L (hence the number of 
modes) constant during the differentiation, in agreement 
with the fact that L = \yA + R?A^ — IJ is constant for a 
mathematically infinitesimal change of R. 

We may obtain a more intrinsic expression for t. With 
{A) = Ap{l+a.), where a is given by eq. ([32]). and iVmodcs = 
^2(2/ + 1), we may rewrite r as 



ksT 

rAp = {A) a r— iVmodos ■ 



(47) 



The quickest way to this result is to keep separate, when 
differentiating with respect to R^, the two terms coming 
from ln{Hi) and \n{l/R^) in eq. The interpretation of 
this equation is not straightforward, because j/cbT is the 
internal energy per mode (not the free-energy per mode). 
Note that the same form for r is also valid in the planar 
case, as shown in [3T] (with further justifications in [TB]). 
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In addition, we have checked that differentiating F 
with respect to Am = 47r[i?(l + (u))]^, which differs from 
47ri?^ by terms oiO{u^), gives a wrong result (in tlie sense 
that it differs from the result obtained by the stress tensor 
method) . 



4.1 Poked vesicle 



In the case where the volume is unconstrained, e.g., a vesi- 
cle poked by a micropipette or by specific proteins, the 
natural choice for the reference sphere (of radius R) is 
the average vesicle's shape. By definition, then, uq^ — 0, 
which implies Am = AttR^. In this case, the Hamiltonian 
in the Gaussian approximation becomes: 



(48) 



where H[ ~ Hi + Ana [ID] . Equation ((55|) , which gives r 
by the stress tensor method, is valid whatever the form of 
Hi. Hence, we need just replace Hi by Hi + 4:K,a, which 
yields: 



I {I + 1) {21 + I) 



(49) 



(To avoid confusions, we denote here by r' the mechanical 
tension in the case of poked vesicles.) Although eqs. (|49p 
and (|36p differ mathematically, it turns out that their dif- 
ference is numerically irrelevant. Indeed the extra term 
4(t/[(Z — 1)(/ + 2)] in the denominator of eq. (|49)) is only 
important for small Vs while the sum is dominated by 
large Vs. 

The free-energy, is given by the same expression as 
eq. with Hi replaced by /f/: 



F' = AwR^a + fceT 



21 + 1 



I 



In (PH'JR^ 



(50) 



It turns out, again, that r' = 9F'/9(47ri?^) exactly. This 
result is satisfying, but at the same time it shows how 
slippery the free-energy approach can be: differentiating 
with respect to the area of the average vesicle is correct 
in the case of poked vesicles but not in the case of closed 
vesicles. The stress tensor method is thus a much safer 
method. 

Our expression for r' differs from that obtained in 
Ref. [50], where the authors considered also a spherical 
membrane without volume constraint. In particular, the 
mechanical tension obtained in that reference cannot take 
negative values. We believe that the discrepancy between 
the two results comes from the omission in Ref. of the 
factors R within the measure. Indeed the factor 1/R^ in 
the logarithm of our eq. ([50)) is absent in the corresponding 
expression (A. 9) of Ref. pO] . 



5 Summary and discussion 

In this paper we have applied the stress-tensor method |16j 
to the calculation of the surface tension r of quasi-spherical 
vesicles, investigating both closed vesicles (volume con- 
straint) and poked vesicles (no volume constraint). In both 
cases we have shown that the more standard method based 
on differentiating the free energy with respect to the pro- 
jected area gives the same results, provided it is performed 
in the right way. 

In order to describe the elasticity of the vesicle we 
have used the simplest Helfrich model, since all its vari- 
ants (ADE. BC, etc.) are equivalent in the quasi-spherical 
limit jl3j . We have made two approximations: i) we have 
replaced the membrane area constraint by a Lagrange 
multiplier approach with effective surface tension cr; ii) we 
have kept all calculations at the Gaussian level. The for- 
mer approximation is justified for small excess areas (see 
sec. 4C of Ref. [13]). To control the validity of the latter 
approximation, we have restricted our investigations to 
relative height variations of the membrane with respect 
to the reference sphere of less than 5%. 

We have calculated the projected stress tensor in spher- 
ical geometry. This tool should be useful in other con- 
texts. We have thermally averaged it in order to calculate 
the actual surface tension r acting tangentially along the 
average membrane surface. Its expression in terms of R 
(vesicle radius), cr, k, k^T, A (cutoff) depends on whether 
the vesicle is closed or poked. However, numerically the 
difference is practically indistinguishable. The flat mem- 
brane limit holds when R > 1 jim. We have re-obtained 
the same expressions for r by the free-energy approach. 
This confirms that the so-called naive measure that we 
have used is valid in the Gaussian approximation (see, 
e.g., Ref. [H]). We have discussed whether the free-energy 
should be differentiated with respect to the area of the 
sphere that has the same volume as the vesicle or with re- 
spect to the sphere corresponding to the average vesicle's 
shape. We have shown that the correct choice depends on 
whether the vesicle is closed or poked, and therefore that 
the stress tensor approach is the safest one. 

It was shown in Ref. [T3] that a (the Lagrange mul- 
tiplier) may become negative. One of our motivation was 
to check whether the actual mechanical tension r could 
become negative or not. This would imply that vesicles 
could sustain inner pressures lower than the outer pres- 
sure (which is impossible for liquid drops). Our analysis 
shows that r may indeed become negative within the va- 
lidity range of our Gaussian analysis. For instance, small 
vesicles with i? = 50 nm can reach negatives tensions of 
— 10~*N/m with an excess area still less than 1% (see 
fig. [3]). The so-called "giant vesicles", with R > l^m, can 
sustain negatives tensions ~ — /cbT^^^/Stt which may be of 
order — 10~^ N/m or — 10~^ N/m depending on the uncer- 
tainty on the value of the cutoff. The possibility to sustain 
negative tensions, or negative pressure differences, might 
be investigated: i) by controlling the outer osmotic pres- 
sure, in the case of small vesicles, or ii) by poking a giant 
vesicle through by a micropipette to which it would adhere 
and gently decreasing its inner pressure. 
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An interesting feature, which might be tested experi- 
mentally by controlling the pressure outside the vesicle, is 
the plateau of fig. [31 One sees that when t reaches a criti- 
cal value Tc < 0, the excess area rises dramatically. As dis- 
cussed above, this behavior (still within the validity range 
of our Gaussian approximation) probably corresponds to 
the transition to oblate shapes. For small vesicles we find 
that it corresponds roughly to TcB? / k w —5 while for 
giant vesicles it is given by Tc ~ —kBTA^/{%'K). Hence, 
applying the Laplace pressure formula AP — 2rc/i?, we 
find that the critical pressure difference yielding the shape 
transition is AP^ « -s\xp[lQK/R^,kYiTA^ /{At:R)]. 

Finally, a possible way to measure the difference be- 
tween r and tr, suggested to us by L. Limozin [35], would 
be to pursue experiments on adhering vesicles similar to 
those of Raedler et al. [31]. In these experiments, the 
tension is deduced both from the fluctuation spectrum 
of the adhering part and from the contact angle "daS via 
the Young-Duprc relation, the adhesion energy being es- 
timated using the distance to the substrate. From our un- 
derstanding, the former method should give a and the lat- 
ter one should give r. By way of illustration, the results 
of Ref. [3T] for a ~ 1.7 x 10~^ J/m^ gave an adhesion en- 
ergy — cos(t9off)] — 1-3 X 10~^ J/m^ significantly larger 
than the value Wa = 8 x 10~^° J/m^ estimated from the 
membrane-substrate separation. The authors of Refs [21] 
and |32j point out that this problem is recurrent. If the 
correct relation is rather t[1 — cos(t?off)] = Wa then we 
obtain r ~ 1.1 x 10~^ J/m^, which is compatible with the 
formula cr — t ~ kBTA^/{STT) [see eq. ([55])] if one takes 
~ 5nm. More experimental studies along these lines 
would be welcome. 



A.I. gratefully acknowledges the hospitality and support of the 
Laboratoire Matiere et Systemes Complexes (MSG) of the Uni- 
versite Paris Diderot-Paris 7 where this work was initiated. 



A Derivation of the stress tensor 

In Sec. 12.21 we explained the general method to derive the 
stress tensor. Here we give the details of this calculation. 
In order to compare cqs. ^ and (|TU)) . one needs to obtain 
59, d4>, 5u, Sug and in terms of 5a. The first three 
can be obtained by identifying the new membrane shape 
u = u + Su with the old one translated: R[l + u{9 + 69 , (fi + 
(50)] e,. (61 + (561,0 + (50) = R[l + u{9 , (p)] Br + Sa. This leads 
to 



torques) ii Sn ■ tg = and 6n ■ = over the border. 
This leads to 



(5a,- ~ R {ug 59 + (50 - 
Sag = Ril + u)69, 
Sa^ = i? (1 + u) sin 9 5(f) 



5u) , 



(51) 
(52) 
(53) 



Let 5n = n{9 + (56*, + 50) — n(6', 0) be the variation of 
the normal, where n is the normal to the shape defined 
by m(0, 0) and n = tg x t^/\tg x with ti = dir. The 
variation of the normal vanishes (implying no work of the 



5ug = \59\{1 + uf + 2ul - (1 + u) 

1 + M L 



ugg\ 



5<j) 



H — : — - (1 + u) U(h cos 6* — (1 + m) ug^ sinl 
sm 9 



(54) 



+ 2u^ug sin ] + (5w We I , 
5u^ = — ^39 [ (1 + li) ii0 cot 9 + 2ugUij, 

- (1 + It) ug^] + (50[ (1 + uf sin^ 9 + 2u^ 

— (1 + u) {ug cos 6* sin 6^ + u^^) ] + (5uu^| . (55) 

These equations, combined with eqs. (|STH55)) . allow us to 
write 5ug and in terms of 5a. Up to order w^, we 
obtain 



69=^-^{l-u + u^), 



(56) 
(57) 



i?sin6l 

5u ~ — [5ar — (1 — u) {Sa^u^ esc 9 + 5agug)] , (58) 
-ft 

5ug = |(5ar (1 - u)ug + 5ag [l + - (1 - u) Ugg] 

+ (5a0 CSC 9 [m0 [(1 - m) cot 6* + ue] — (1 — u) ug^] | , 

(59) 

6u(j, = — |(5ae [u^ [(1 - u) cot 9 + ug] - (l - u) ug^ 
+ 5ar (1 — u) W0 + 5a^ esc 9 [ sin^ + 
— (1 — u) (u00 + Me cos6'sin0) ] I . (60) 

These expressions are to be inserted into eq. ©. Note that 
it is necessary to expand h at 0{u^) in order to obtain 
dh/dui and dh/duij consistently at 0{u^) in eq. This 
means adding 

/13 = — K sin {4 esc"' 9 U(f, ug ug^ + 2ul ugg 
— 2u0 csc^ 9 esc^ 9 — ug cot 9) 
+ u [ugg +ug{2 + cot^ 9) + 2u| csc^ 9 + it|^ esc"' 9 

+ 2ug Ugg COt 9 + 2ug U^^ CSC^ 9 COt 9 + 2ugg CSC^ 6*] 
+ 2u^ [ugg + Ug COt 9 + U^^ CSC^ 9) } , (61) 

to eq. ^ before calculating the derivatives. Comparing 
eqs. ^ and pH)) . we obtain finally 



see eqs. 

These expressions arc valid up to 0{u'^). 
A.l Planar limit 

In order to check the expressions (|^ - |57|) , we consider the 
limit i? — >■ oo, which should yield the stress tensor for a 
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See = ^ { -^^'^ 6 (2 + 2u + u1 csc^ — Ug) + k [u1^ csc^ - uee sin + 2ue {ue,p,p esc + ueee sin 6)] 
+^080 6* [2it0 — ug cos^ 6* — 2it0</, (1 + ug cot + 4it (it^,/, — itee sin^ 

+2k [uee (sin 6 + ug cos ^) + (2u — l)ue cos 5] } , (62) 
S4,e = —Raugu^, + {2ug^ — cot 9) + ^ (4u cot 6 — u^^ cot S csc^ 6^ + iiew^ — 4tt ue^ — u^^ ug^ csc^ 6 

— We ue^ cot 61 + it^ ite^^ csc^ 6 + 2u^ uee cot 6 — ue4,uge + u^ueee) , (63) 
Srg = Raug sin — [(1 ~ 2it) ite (2 sin 6 — esc 0) + it^e cos 6 — 2u^^ cot 6* esc 9 + ue,^^ esc 6* + ugee sin ^] 

— — [2?i| cot 6^ CSC 8 — ue cos S + 4u u^^ cot esc 6 — 2u uge cos 6' 

— 2u ue00 CSC 6 — 2u ueee sin ^ + u^^ ue esc 9 — 3itg itge sin ^] , (64) 

lr2/ 22 2\ r24: 2 2/ 2 \"l 

£'<;)0 = ^ |-R CT (^2 + 2?i — it0 esc 6 + ugj + k y—u^^ esc 6 + ugg + 2u^ CSC 9 (^11^0^ CSC 6 + ues^jj 

-fftcsc^ 9 [2u0 + Ug (3 sin^ ^ — l) + 2u<^0 (1 — ug cot 6* — 2it)] + 2k (1 — 2ii) tie cot ^ — 2k (1 — 2u) uee 

-|-2ku0 ue0 csc^ 6^ cot 6'} , (65) 

Sg^ — —Raug 110 esc 9 H ^ — [ {u^ cot 6 — ug^,) (—2 + 4ii + itee + u,f,^ csc^ 6*) + ue u^ (l + csc^ 6*) 

+ ite(ti000 csc^ 61 + itee^)] , (66) 
Srcj> = [^-^i (-R^cr — 2k + 4Kit + 3it00 esc^ 9 + 3ue cot 6^ — uee) — k (1 — 2u) (ue^ cot 9 + ugg^ + u^^^ csc^ 6*)] . (67) 



flat membrane [TU]. We consider a general point (0, 0) on 
the sphere and we define a local system of cartesian coor- 
dinates (x, y, z) such that da; = Rd9^ dy = i?sin0 and 
dz = dr. We want to determine the stress tensor S' in 
these new coordinates. The component of the elementary 
force dfa along a general direction a € {x,y,z} exerted 
through a cut perpendicular to x reads dfa = S'^x di ~ 
Sa0 d(f>, with the correspondence d£ = R sin 9 d(j). Thus 
we have S'^x = Sae/{Rsin0). Likewise, E'^y = Satj,/R 
since in this case d£ = Rd9. In the limit of large R, 
the plane tangent to the sphere at the point ((/), 6) be- 
comes the reference plane of the membrane, and the height 
of the membrane over this plane is h{x,y) — Ru{9,(j}). 
Hence, ug = Ru^ = h^, = R sin 9 Uy = sin 9 hy, 

Uee = R^ Uxx = R hxx, u^<p = R^ sin^ 9 Uyy = R sin^ 9 hyy, 
etc. Keeping the terms that are dominant in the limit 
i? — )■ (X) , we obtain 



y — 

^xx 



nhxdx^^h , 



yx — —dhxhy — nhxyS/ h + Khydx^ h, 



^zx = (^hx - ndxV^h, 



Eyy = a 



{hi -hi) 



"'yyJ 



+ nhydyV^h , 
^xy = —<yhxhy — nhxyV^h + Khxdy\7^h, 



Sry = crhy — KdyS/^h, 



(68) 
(69) 
(70) 

(71) 
(72) 
(73) 



which agree with [TU] . 



B Correlation functions and relations among 
them 

There are five fundamental correlation functions, which 
we have calculated explicitly: 



see eqs. 



The other correlation functions either vanish or may be 
deduced from them, as explained below. Using the rela- 
tion: 



gp-iy/ 



i,n+l 



Qpyn 



d(j)P 



(74) 



which holds also in the presence of derivatives with re- 
spect to 9, one deduces the following rule: when averaging 
a product of two terms one may pass a derivative with 
respect to cj) from one term to the other one while multi- 
plying by —1. hence. 



(u^Ugg) 



(uu^g) , 
{ug U^e) 



(75) 
(76) 
(77) 



An interesting consequence is that averages implying an 
odd number of derivatives with respect to (p vanish: 



{uu^) = , 
{ug u^) = , 
(u0 u^^) = , 



(78) 
(79) 
(80) 
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{ue ue^) = , 
{ubo Ug^) = , 
{Ucficl, U0^) = . 



(81) 
(82) 
(83) 



Starting from the addition theorem for spherical harmon- 
ics, one may show by suitable differentiations that the sum 



m— — / 



Qpyr. 



del' 



OOP 



vanishes for k + p odd, implying: 

{u Ue) = , 
{uugee) = . 



(89) 



(90) 
(91) 



Note that this does not hold when derivations with respect 
to (j) are also involved. 

Starting again from the addition theorem, one may 
show that 



m— — l 



dQk+l 



m— — l 



86'' 



d0P 



(92) 



It follows that, when averaging a product of two terms one 
may pass a derivative with respect to 6 from one term to 
the other one while multiplying by —1. This holds only, 
however, in the absence of derivatives with respect to 0. 
As a consequence. 



{U0 Ue) 

{ue Uggg) 



{uuee) , 
{uee uee) 



(93) 
(94) 



Finally, one may also use the fact that AY[^{9, (p) = 0, 
where A is the laplacian in spherical coordinates, to obtain 

{ue uee) = cot 6{u uee) - csc^ 8{ue u^^) . (95) 

Since {ue uee) vanishes, on obtains 



[Ue 



= — sin cos 9 I 



(96) 



In conclusion, eqs. ([MH881) . together with eqs. 
eqs. (|90H9T|) and eq. ((96|) give all the correlation functions 
of the derivatives of u. 
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